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Abstract. In this paper we present a non-Gaussian integral based on a cubic polynomial, 
instead of a quadratic, and give a fundamental formula in terms of its discriminant. It 
I gives a mathematical reinforcement to the recent result by Morozov and Shakirov. We also 

' present some related results. This is simply one modest step to go beyond the Gaussian but 

it already reveals many obstacles related with the big challenge of going further beyond the 
Gaussian. 
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^ ■ 1 Introduction 

The Gaussian is an abbreviation of all subjects related to the Gauss function e~^^^ +qx+r) j-j^g 
the Gaussian beam, Gaussian process, Gaussian noise, etc. It plays a fundamental role in math- 
^ . ematics, statistics, physics and related disciplines. It is generally conceived that any attempts to 

^ I generalize the Gaussian results would meet formidable difficulties. Hoping to overcome this high 

wall of difficulties of going beyond the Gaussian in the near future, a first step was introduced 
in [l]. This paper is its polished version. 
CN ■ In the paper [2j the following "formula" is reported: 

(N ■ 

where D is the discriminant of the cubic equation 

ax^ -|- hx^ -|- cx -|- (i = 0, 
and it is given by 

B = \?(? + \Mhcd - W - M?d - Tla^d^. (2) 

The formula dH) is of course non-Gaussian. However, if we consider it in the framework of the 
real category then ([1]) is not correct because the left hand side diverges. In this paper we treat 
only the real category, and so a, 6, c, d, x, y are real numbers. 

Formally, by performing the change of variable x = tp, y = p for ([1]) we have 

l.h.s. oi& = JJ e-P'i''''+''"+^'+^)\p\dtdp = y {/ e-i''''+'""+'^+''y\p\dp^ dt 



^3 



Idle da I - — , — =dt 

^{at^ + ht'^ + ct + d)\^{at^ + bt"^ + at + d) 



by the change of variable a = -\/ at^ + ht"^ + ct + dp. 
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The divergence comes from 

3 

\cr\e~'^ da, 
while the main part is 

1 

zdx 



I 



{ax^ + 6x2 _|_ _|_ ) j ^ (ax3 + bx"^ + cx + d) 
under the change t — )■ x. As a kind of renormahzation the integral may be defined like 

J Jr2 Jr I -y/ {ax^ + + cx + d)\'^ (ax^ + 6x2 + cx + d) 

However, the right hand side lacks proper symmetry. If we set 

F{a, b,c,d)= [[ e-('^^'+^^'^+'^^^'+'^3'')dxciy, 

where Dfi = [—R,R] x [—R,R], then it is easy to see 

F{—a, —6, — c, —d) = F{a, 6, c, d). 

Namely, F is invariant under Z2-action. This symmetry is important and must be kept even 
in the renormahzation process. The right hand side in the "definition" above is clearly not 
invariant. Therefore, by modifying it slightly we reach the renormalized integral 

Definition 1. 

t [ [ e-("^'+^^'s'+'^^2''+'^2^')dx(iy t = [ ^ (3) 

J Jir2 Jjj ^ (a2;3 _)_ 5a;2 -\- cx + (i)2 

We believe that the definition is not so bad (see the Section [3D. 

In the paper we calculate the right hand side of ([3]) directly, which will give some interesting 
results and a new perspective. The result gives a mathematical reinforcement to the result [2] 
by Morozov and Shakirov. 

2 Main result 

Before stating the result let us make some preparations. The Gamma function T{p) is defined 
by 

T{p)= / e-^x^-^dx ip>0) (4) 
Jo 

and the Beta function B{p, q) is 

B{p,q)=[ xP'\l-xy^^dx {p,q>0). 
Jo 

Note that the Beta function is rewritten as 



oo x^-i 



B(p,q) = / -; ——dx. 

' (l + x)P+'' 







See [3] for more detail. Now we are in a position to state the result. 
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Fundamental formula. 

(I) For D < 



where 



' 2' 6 



(II) For D>0 



I 



^ d^=^^ (6) 



^{ax'^ + 5x2 + CX + VD 

where 



3' 3^ 

(III) C_ and C+ are related by C+ = -v/SC^ through the identity 

Our result shows that the integral depends on the sign of D, and so our question is as follows. 

Problem. Can the result be derived from the method developed in ^? 

A comment is in order. If we treat the Gaussian case (^g-i°-^^+^^v+'^y^)'j then the integral is 
reduced to 

1 , 27r 

— K ; dx = . (8) 

ax2 + bx + c 

if a > and D = b'^ - 4ac < 0. Noting 

^_V^V^ r(i)r(i) ^/ii 



1 r(i) \2' 2 

should be read as 



ax2 + 6x + c \/—D 

3 Proof of the formula 

The proof is delicate. In order to prevent possible misunderstanding we present a detailed proof 
in this section. 

Proof of (I). We prove dS]) in case of D < 0. 

First we consider the special case where a = in the cubic equation ax^ + bx^ + cx + d. 
Namely, we calculate the integral 



/ 

Jm 



{/(5x2 + cx + (i)2 



dx. 
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Noting -D = b'^{4bd - c^) > we obtain 



/ 



zdx 



zdx 



zdx 



1 



c_\2 I _ c2 



' b a? 



1 



62 A ^(a;2 + r2)2 



dx 



z=rj/ 1 



zdx 



T 



1 



^Pryo ^(y2 + 1)2 

1 

1 



F A v^^(y2 + 1)2 
1 



62ryo ^(x + l)2 2Vi 



62^70 (X + I)i vW2 y62(45^_c2) 



(9) 



Now we consider the general case of a 7^ 0. From the condition Z? < there is (only) one real 
root of the cubic equation ax'^ + 6a;2 + cx + d = 0. Let us denote it by a. Prom the equation 

ax^ + bx^ + cx + d = {x — a)(ax^ + kx + I), aa^ + ba^ + ca + d = 

we have easily 

b = k — aa, c = I — ka, d = —la. (10) 
First we assume a = 0. In this case d = and 
ax^ + bx^ + cx + d = x(ax^ + bx + c) . 

Then 



zdx 



Y^x2(ax2 + 6x + c)2 Jo Y^x2(ax2 + 6x + c) 



1 



.1 rO 



2 V 



41 + 



1 



Y^x2(ax2 + bx + c)2 



=dy + 



i/{cy^ + by + af J-00 ^ {cy^ + by + a) 



y \y y 
dy 



fdy 



{/(cy2 + 6y + a)2 
Next, let us calculate the case a 7^ 0: 



-y/(x — a)2(ax2 + kx + 



dx 



x=y+a 



dy 



^y^{a{y + af + k{y + a) + lY 
dy 



^y'^iay'^ + {2aa + fc)y + {aa^ + ka + l)Y 
^J{aa'^ + ka + /)2{4a(aa2 + ka + l) - {2aa + k)"^} 



^{aa^ + ka + lY{Aal -k^y 



(11) 
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Key Lemma. From (jlOp the following equation holds 

(oQ^ + ka + l)'^{4al - k^) = 21c? + 4ac^ - XMhcd - }?(? + 46^d = -D. (12) 

The proof is straightforward but tedious. 

Therefore, from both (HID and (1121) we obtain the formula 



y[x-Oif(ax'^^kx^lf 

Proof of (II). We prove ([6]) in case of Z) > 0. Let us start with the evaluation of the following 
integral 

zdx 



/m {/x2(x — a)2 
for Q > 0. Then 

fOO 



[ — ==dx = I , „ . .„ dx + / „ , „ , ,„ dx 

Jr \/x'^(x ■ ' " 



1 

-y/x2(x — a)2 j-oo y^x'^{x — a)2 Jo Y^a;2(2; — a)^ 



— ^dx + / — ^dx, (13) 

-y/x2(x + a)2 Jo ^^^(x — a)^ 

where the change of variable x — )• — x for the first term of the right hand side was made. 
Each term can be evaluated elementarily: 

/ — =ax = — ^ / — =at = 03/ = a - - 

Jo i/x^ix + ay ^Jo ^W+lF ^0 (i + l)l V3 3 

while 

^ ^ poo 2 

^x2(x-a)2 io syt2(t- 1)2 

a"5 J / , ^ dt+ [ , ^ dtl 



1 



a~3 < / =dt + / zdt 

/o {/t2(l-t)2 ^t2(i _ 1)2 



1 



2a" 



3 / - d.t = 9.^-3 / t~3(l-t)-3dt 



^t\l-t) 
1 1 
3'3.' ' 



where we have used 



h ytHt - 1)2 7i ./ 1^ (1-.)' V sV 7„ y»2(i-s)2 io yt\i-ty 

From (1131) we have 



"71 c'2 



f 1 , 1 A 35(i i) 

/ :dx = Sa^^B -, - = 

Jr {/x2(x-q)2 V3 3; ^ 
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Now we consider the special case a = in the cubic equation ax^ + hx^ + cx + d. Then by 
D = h'^{c^ - Ahd) > we obtain 



1 



+ CX + d)' 



a^ = ^>0 1 



zdx 



1 



OO 3 



zdx 



[h{x + 



_c_\2 _ c^—Ahd \ ' 
2b > 4b j 



1 



zdx 



46 



=x+a 1 



^(x2 - a2)2 



zdx 



\/{x — aY{x + a)' 



zdx 



IP 1 3i?(i,i) 



/62 ^ 



3i?(i,l) _ 3ij(i,|) _ 
\/2a^ v^i^ V62(c2 _4ftrf) 



(14) 



Next we consider the remaining general case of a ^ 0. From the condition D > there 
are three real solutions in the equation ax^ + bx'^ + cx + d = 0. We denote one of them by a. 
Remember the relations b = k — aa, c = I — ka, d = —la from the equation 

ax^ + bx^ + cx + d = {x — a) [ax^ + /cx + Z) , aa^ + bo?' + ca + d = 0. 

First we assume a = 0. Then 

ax'^ + 6x^ + cx + d = x(ax^ + bx + c) 

and from D = cP'ilP' — Aac) we have 



zdx 



■^7x2(0x2 + 6x + c)2 Jo {/x2(ax2 + 6x + c) 



ifix + 



-1 rO 

y 



OO 3 1 



y2 



+ 



■^7x2(0x2 + 6x + c)2 

(4/ 

2 



dx 



y \y y 



y 



fdy 



^{cy'^ + by + a)2 

ITU (c^-fe; fe^-c; a^-d) 35 (|, |) 



Vc2(62 - 4ac) 

For the case a 7^ we obtain the formula 



Ay(cy2 + by + a)2 



dy 



^(cy2 + by + a)2 



(15) 



/ 

■718 



^(x — a)2(ax2 + kx + /)2 



dx 



1 



GUI 



{/y'^{ay'^ + {2aa + A;)y + (aa2 + A;a + /)}2 

1 1 

,3' 3 



{/y2{a(y + q)2 + k{y + a) + 0^ 
dy 



dy 



3i?(i 



y(aa2 + to + 0^{(2aa + k^ - Aa{aa'^ + ka + I)} 

3i3(i,i) m^Bjll) 



^{aa^ + ka + lf{k'^ - Aal) 



Proof of (III). We prove the relation ([7]). Let us make some preparations. For the Gamma 
function ([4]) there are well-known formulas (see for example [3]) 



B{x,y) 



r(x)r(y) 

r(x + y) 



(x,y > 0), 



(16) 
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T{x)T{l-x) = . r N (0<x<l), (17) 
sin(7rx) 



psp is called the Legendre's relation. In the formula we set x = 2/3, then 
Multiplying both sides by r(l/6) gives 



3/ \2J \6J \3 



r(if r(i) 



r(i)r(i)r(i) Vs'e 
T^r I 12' 6 



sm 



where we have used formulas (jl6|) and (jl7|) several times. 
The proof of ([7]) is now complete. 



4 Renormalized integral revisited 

In this section let us check whether the renormalized integral ([3]) is reasonable or not by making 
use of the results in the Section [3l 

In the introduction we introduced the following integral defined on Dji = [—R, R] x [—R, R] 



F{a, b,c,d)= II e-('''''+'''''y+'^y'+''y')dxdy. 
J J Da 



For this it is easy to see 



Dr 

d d 9 d\ 
dbdb-d^d-c) ^(«'^'^'^) 

= // {x^y ■ x'y - x' ■ xy2)e-K+''-'?'W+d?/^)d^^y = q, 

J JDr 
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.d-cd-c- dbdd^ i^(a,6,c,d) 



Dr 



(19) 



On the other hand, if we set 



F{a, b, c, d) 



zdx 



C± 



i/{ax^ + 6x2 + CX + df \/±D' 
then we can also verify the same relations: 



d d d d . ^ 



d d d d \ , „ 

56 56-9^a^j-^('^'^'^'^) = °' 



\dcdc dhdd)-^^"'^' 



c,d) = 0. 



(20) 



Verification by hand is rather tough, but it can be done easily by use of MATHEMATICA0. 
Prom ()19p and (j20p we can conclude that our renormalized integral ([3]) is reasonable enough. 

5 Discriminant 

In this section we make some comments on the discriminant ([2]). See [1] for more details ([1] is 
strongly recommended). 
For the equations 



f{x) = ax^ + hx^ + cx + d, f'{x) = 3ax^ + 26x + c 
the resultant R{f, /') of / and /' is given by 



(21) 



a h c d Q 

a b c d 

3a 2b c 

3a 26 c 

3a 26 c 



RifJ') 



It is easy to calculate (j22p and the result becomes 



1 



(22) 



-i?(/, /') = 27a^d^ + Aac^ - ISabcd - b^c^ + Ab^d = -D. 
a 

On the other hand, if a, /3, 7 are three solutions of /(x) = in (j2ip . then the following 
relations are well-known 



6 
a 



a + /3 + 7 
From these it is easy to see 

a + /3 + 7 



a(3 + a7 + /37 = — , 
a 



6 
a 



a^ + P^ + 7^ 



2ac 



d 
a 



a^ + j3^+ 73 



6^ + Sa^d - 3a6c 



4 ^4 4 6^ + 4a26d + 2o2c2 - 4a62c 
+ + 7^ = 2 . 



^The author owes the calculation to Hiroshi Oike. 
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If we set 

A = (a-/3)(a-7)(/3-7) 
the discriminant D is given by 



D = a'A^. 



Let us calculate directly. For the Vandermonde matrix 



V 




\V\ 



-A 



we obtain by some manipulations of determinant 



A^ 



(-\V\f = \V\\V^ 



\VV^ 

b 
a 

- 2ac 



3 a + /3 + 7 + /32 + ^2 

a + /? + 7 + ^2 _^ ^2 Q,3 _^ ^3 _^ ^3 

a2 + /32 + 72 a3 + ^3^^3 „4 ^ ^4 ^ ^4 

6^ - 2ac 

^2 

+ 3a'^d- 3abc 



a 

- 2ac 



2b 



b^ + 3a^d - 3abc b^ + Aa'^bd + 2a^c^ - iab'^c 



2ac 



a 
2c 
a 



Sad — be 



26^ + 2ac Sad - 36c Aabd + 2ac^ - 2b'^c 



3 


b 


6^ - 2ac 




3 


b 


2c 




a 


a2 






a 


a 


2b 


2c 


3a(i — 6c 




26 


2c 


Sad + 6c 


a 


a 


a2 




a 


a 


a2 


2c 


Sad — be 


abd + 2ac^ - 62c 




2c 


Sad — be 


-2bd + 2c^ 


a 


a2 


a3 




a 


a2 


a2 



3 

2 



6 2c 
a a 
6^ — 3ac 6c — 9ad 



3a2 



6c — 9ad c 
— — — 2 



3a2 

2 - 36(i 



3a2 



3a2 



6^ — 3ac 6c — 9ad 



3a2 
6c — 9ad 



3o2 
-2 - 36d 



3a2 



3a2 



1 -1 



= — — {(6c - 9adf - 4(6^ - 3ac)(c2 - 36d)} . 

Qi o 

This result is very suggestive. In fact, from the cubic equation 



ax'^ + bx +CX + d = 



we have three data 

^ = 62 - 3ac, B = bc- 9ad, C = - Sbd, 
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and so if we consider the quadratic equation 

AX"^ + BX + C = 

then the discriminant is just — AAC. This is very interesting. 
Problem. Clarify the above connection. 
As a result we have 

D = — {{bc- 9adf - 4(6^ - 3ac){c^ - 3bd)} = b'^c^ + ISabcd - Aac^ - Ab^d - 27a^d^. 
3 

6 Some calculations 

In this section we calculate some quantities coming from the integral. 
The expectation value (x^) is formally given by 

(x^) = = _^ log J / / e<''-'+^-'y+c-y'+dy')dxdy \ , 



SO renormalized expectation values (x^)rn5 {x'^y)RN, {xy'^)RN, {y^)KN are defined as 
Definition 2. 



{x')kn = log {t /^^ e-^'^^'+'^'y+^^y'+'^y'^dxdy t} , 

{x'y)RN = -^log |t 11^^ e-iax'+bx^y+cxy^+dy')^^^y j| 
{xy'h^ = log |t 11^^ e-^ax^+b^'y+cxy-+dy^)^^^y j| 

(y')RN = log It // J 



dd 

From the integral forms ^ and ^ it is easy to calculate the above. Namely, we have 

, ISbcd - 4c3 - 54a(i2 , ^ ^ 2bc^ + 18acd - 12b^d 
{x )rn = ^ , {x y)RN = ^ , 

, 26^0 + ISabd - Uac^ , n, 18a5c - 46^ - Ma^d 
{xy )rn = ^ , {y )rn = ^ , 

where D = b'^c^ + ISabcd - Aac-^ - Ab^d - Tla^d^. 

We can calculate other quantities like {x^y)Ym or (x^y^)RN by use of these ones, which will 
be left to readers. 

7 Concluding remarks 

In this paper we calculated the non-Gaussian integral ([T]) in a direct manner and, moreover, 
calculated some renormalized expectation values. It is not clear at the present time whether 
these results are useful enough or not. It would be desirable to accumulate many supporting 
evidences. Some application(s) will be reported elsewhere [5]. 
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At this stage we can consider a further generahzation. Namely, for the general degree n 
polynomial 

f{x) = aox" + aix""^ H h a„_ix + an 

the (non-Gaussian) integral becomes 

1 



zdx. 



(23) 



The discriminant D of the equation f{x) = is given by the resultant R{f,f') of / and /' 



like 



1 . n(n-l) 

ao 



where 



D = {-l)'^R{fJ-')/ao 



R{fJ') 



ao ai 
ao 



nao {n — l)ai 
nao 



(n — l)ai 



On— 1 

ao ai 

O-n-l 



a-n-i 



nao ~ l)fli 



see pT]) and 

For example, for n = 4 



R{f,f') 



ao 


ai 


02 


as 


04 











ao 


ai 




as 


04 











ao 


ai 


a2 


as 


04 


4ao 


3ai 


2a2 


as 














4ao 


3ai 


2a2 


as 














4ao 


3ai 


2a2 


as 














4ao 


3ai 


2a2 


as 



and 



QQ QQ 0/1 /10 OOO /I 

= 2560904 — 40^03 — 270903 — 270^^04 — 128090204 + 0^^0203 + I600O2O4 

oQQo 00 00 00 

— 4O0O2O3 — 4O1O2O4 + I44O0O2O3O4 — 6ooO]^0304 + 144ooO;^0204 — 1920oOiOsa4 



+ 1800010203 + 180^020304 — 80ooaio|o304, 



and for n = 5 



RifJ') 



ao 


01 


a2 


as 


04 


05 














Oo 


01 


a2 


as 


04 


05 














Oo 


01 


a2 


as 


04 


05 














Oo 


01 


a2 


as 


04 


05 


5oo 


4oi 


302 


203 


04 

















5oo 


4oi 


302 


203 


04 

















5oo 


4oi 


302 


203 


04 

















5oo 


4oi 


302 


203 


04 

















5oo 


4oi 


302 


203 


04 
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and 

D = SUba^al - 2500a|]aia4ai - 3750a|]a2a3a5 + 2000a|]a2ala^ + 2250a|]aia4a^ 

— 160000030405 + 256aj]o4 + 2000aoof 030^ - 50000^0405 + 2250ooOi02o| 

— 2050000102030405 + 1600001020405 — 9OOO0O1O3O5 + 10200001030405 

— 192aoOi03a| — 90000020405 + 82500020305 + 5600002030405 — 128ooa2a| 

— 6300002030405 + I44O0O2O3O4 + IO8O0O3O5 — 27O0O3O4 — 1600ooO]^0205 

00 'i'i 

+ 1600901030405— 8600010405 + 10200001020405+ 560000^020305— 74600O1O203O4O5 

OA O O O 'X 

+ 144000^^0204 + 24000^^030405 — 6000^^0304 — 63O00O1O2O3O5 + 2400O1O2O4O5 

'\ A "X ^ 

+ 356000102030405 — 800001020304 — 720001020305 + 180001020304 + 108000205 

A. A '\ "X O ^"-i AO 

— 7200O2O30405 + I600O2O4 + 1600020305 — 400020304 + 2560105 — 19201020405 

A AO A A O "X O "X 

— 12801O3O5 + 14401030405 — 270104 + 14401020305 — 601020405 — 800102030405 
+ 180^020304 + 160^0305 — 40^0304 — 2702^0205 + 180^^02030405 — 40^0204 

^223 , 2222 

— 4O1O2O3O5 + 01020304. 

However, to write down the general case explicitly is not easy (almost impossible). 

Problem. Calculate (|23p for n = 4 (and n = 5) directly. 

The wall called Gaussian is very high and not easy to overcome, and therefore hard work will 
be needed. 

Recently the subsequent paper [6] by Morozov and Shakirov appeared. Our works are deeply 
related to so-called non-linear algebras, so we will make some comments on this point in a near 
future. As a general introduction to them see for example [7]. 

The author thanks referees and Hiroshi Oike, Ryu Sasaki for many useful suggestions and 
comments. 
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